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SIMMARY 


A small total-pressure tube resting against a flat-plate surface 
was used as a Stanton tube and calibrated as a skin-friction meter at 
various subsonic and supersonic speeds. Laminar flow was maintained 
for the supersonic runs at a Mach number M,,,, of 2. At speeds between 
Moo = 1.33 and- M» = I. 87 , the calibrations were carried out in a turbu- 
lent boundary layer. The subsonic flows were found to be in transition. 

The skin- friction readings of a floating-element type of balance 
served as the reference values against which the Stanton tube was 
calibrated. 

A theoretical model was developed which, for moderate values of 
the shear parameter t, accurately predicts the performance of the 
Stanton tube in subsonic and supersonic flows. A "shear correction fac- 
tor" was found to explain the deviations from the basic model when r 
became too large. Compressibility effects were important only in the 
case of turbulent supersonic flows, and they did not alter the form of 
the calibration curve . 

The test Reynolds numbers, based on the distance from the leading 
edge and free-stream conditions, ranged from 70,000 to 875,000. The 
turbulent -boundary-layer Reynolds numbers, based on momentum thickness, 
varied between 65O and 2,300. Both laminar and turbulent velocity pro- 
files were taken and the effect of pressure gradient on the calibration 
was investigated. 


INTRODUCTION 


In recent years there has been a renewed interest in the study of 
boundary- layer separation, especially at high speeds. Measurements of 
local shear forces are important in such investigations, but absolute 
skin-friction measuring devices of the floating-element type of balance 
are not suitable; the pressure gradients in the flow give rise to forces 
across the face of the floating element which may be of the same order 
of magnitude as the shear stresses to be measured. A device of much 
smaller characteristic length is needed to overcome this difficulty. 
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Such an instrument is the Stanton tube I half -pitot surface tube). 

It consists essentially of an obstruction wiuh a pressure orifice on 
its front face, the height of which is small compared with the boundary- 
layer thickness. In particular, a very small total-pressure tube resting 
against the surface can serve as a Stanton tube. The ability of the 
Stanton tube to tolerate fairly steep pressure gradients is verified by 
recent studies of shock-wave— boundary- layer interaction (ref. 1) in which 
shear distributions were measured through the interaction region with 
the tube whose calibration is reported herein. 


Until this investigation no reliable calibration of the Stanton 
tube had been made at supersonic speeds. Some preliminary investigations 
were conducted in turbulent supersonic flow \,ref. 2) and they yield d 
inconclusive results. 

The Stanton tube was invented in 1920 and used by Stanton at the 
National Physical Laboratory in Teddington, England. He built a half- 
pitot total-pressure tube for his studies of boundary layers in turbu- 
lent pipe flow (ref. 3). Fage and Falkner (ref. 4) used the Stanton tube 
in 1950 to measure skin friction on airfoil shapes. Taylor (re . 5J 
summarized the existing results and performei additional experiments 
in 1938, essentially under Stokes flow conditions. He formulated the 
calibration problem in terms of two parameters and showed that all 
available data were in reasonable agreement. Further experiments, 
involving the use of the Stanton tube were trade by Fage and Sargent 
(ref. 6) and Cope (ref. 7 ). In both investigations an attempt was made 
to measure supersonic skin friction using sibsonic calibration data. 

The validity of this extension of the calibration is doubtful. 

A significant advance in the calibraticn of the Stanton tube is 
made possible by the availability of an absclute shear meter which 
replaces computed skin friction by measured values. Direct surface 
shear measurements were first taken in 1929 by Kempf (ref. 8) who used 
a floating element on the bottom of a barge to obtain skin- friction 
values at high Reynolds numbers. Not until 1951, however, were there 
any skin-friction meters suitable for high-speed work. Dhawan (ref. 9) 
in 1951 developed a floating- element balance which could be used at 
high speeds. Hakkinen (ref. 2) used an improved version of Dwahan 
instrument to obtain reference shear values and calibrated a Stanton 
tube in high- subsonic and low-supersonic flow. His results for subsonic 
speeds are in reasonable agreement with the older data. The supersonic 
measurements, however, show appreciable dev: .at ions and suggest compres- 
sibility effects. In 1954, on the basis of extensive experiments which 
suggested a unique pressure- shear relation, Preston (ref. 10) proposed 
the use of round pitot tubes resting againsu a flat surface as skm- 
friction meters. Stevenson (ref. ll)in 1955 obtained results which are 
in substantial agreement with Preston's data. 



5 


The object of the present investigation is to extend the Mach num- 
ber range of the Stanton tube calibration, to determine when compres- 
sibility effects become important, and to propose a theoretical model 
to explain the calibration in various flow regimes. 

This investigation was carried out at the Massachusetts Institute 
of Technology under the sponsorship and with the financial assistance 
of the National Advisory Committee for Aeronautics. The authors wish 
to express their thanks to the following persons: Dr. Jeffrey Tinkler 

who designed the basic floating element, Dr. Marten Landahl who graciously 
helped in making the element a working device, and Messrs. Kenneth 
Golden and Isaac Greber who performed the supersonic turbulent experiments . 


SYMBOI£> 


A proportionality constant, defined in equation (2l) 

a speed of sound 

B proportionality constant, defined in equation (22) 

Cj. local skin-friction coefficient 

h characteristic height of the Stanton tube 

M Mach number 

P pressure 

P = p/h 2 

P T Stanton tube pressure 



r Reynolds number based on distance from leading edge and free 

stream conditions 

R 0 Reynolds number based on momentum-thickness and free-stream 

conditions 

s density, nondimens ionali zed through division by wall density 
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T = 

U 

U* 

u = 
u' 

V 

V = 

v' 

X 

X = 

V 

y = 

a 

6 

7 

e 

n 

0 

P 

p 


t ■/hP’ 

velocity in X-direction 

friction velocity in laminar sublayer. 



U/( T v»K) 

- v lfQ^ 

velocity in Y-direction 

v /( T » h /iv) 

= v /f?V 

physical coordinate in streamwise direction 

x/h 

physical coordinate normal to X 

Y/h 

ratio of laminar sublayer thickness to probe height, e/h 
boundary -layer thickness 
specific-heat ratio 
thickness of laminar sublayer 
normal coordinate defined in equation 

momentum thickness of boundary layer. 


(12a) 



coefficient of viscosity 
density 



T E 


shear parameter 
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shear stress at wall. 



T e = U X e K 

Subscripts: 

oo refers to free-stream conditions 

o refers to static conditions 

c corrected 

h refers to conditions at upper edge of Stanton tube 

v refers to conditions at vail 

t refers to conditions that would have existed at upper edge 

of Stanton tube if velocity profile remained linear up to 
Y = h 


EXPERIMENTAL APPARATUS 


Wind Tunnels 

The experiments described herein were conducted in two wind tunnels. 
The subsonic and laminar supersonic runs were made in an 8- by 8- inch, 
continuous -flow supersonic tunnel in the M.I.T. Gas Turbine Laboratory. 
This is a vertical flow tunnel with a fixed nozzle of nominal Mach num- 
ber 2 in which the air flows upwards. 

The turbulent supersonic tests were made in a vertical wind tunnel 
having a 2^ by 2^ - inch closed circuit with variable Mach number and 

downward flow; this tunnel is also located in the Gas Turbine Laboratory. 
The Mach number is varied by means of a sliding nozzle block (ref. 12). 
Wall boundary -layer removal ahead of the test plate is effected by a 
suction duct connected to a steam ejector. 

The pressure level in both tunnels can be changed by means of steam 
ejectors. The stagnation pressure, which controls the test Reynolds 
number per unit length, can be varied from 0.2 to 2 atmospheres. 

The two tunnels are shown in figures 1 and 2. 
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Floating Element 


The floating-element skin-friction meter is shown in figure 3. 

Two configurations, one for each wind tunne.., were used in the present 
investigation. Both differed only slightly from the basic model used 
bv Hakkinen (ref. 13). In a horizontal tunnel, movement of the floating 
element (fig. 3(a), part ®) is due to the shear forces only in the 
absence of chordwise pressure gradient. In a vertical tu™ 1 © accoun 
must be taken of gravity. A counterweight (§) glued to the balance 
arm (C) serves to achieve a "neutral zero piint. Another departure 
from the basic model is the modification of the dashpot pamper (DJ . 
Originally the hollow cylinder © containei the damping fluid and a 
plunger which served both as a counterweight and as the dashpot piston. 

It was found, however, that surface tension and capillary suction due 
to the small clearance between the plunger and the cylinder walls pre 
vented, the use of that configuration. 

In addition, the two modified configur it ions of the floating element 
differed from each other in two respects: 


(l) The vertical air flow in the two tunnels was in opposite direc- 
tions, so that whenever the device was transferred from one test sec- 
tion to the other it had to be inverted. 


(2) Since the flow in the variable Mach number tunnel was turbu- 
lent, stiffen springs © (fig. 3(a)) replaced the ones used in the 
larger tunnel in order to enable the instrument to sustain bigger loads. 


It should be noted that the whole device was designed to be as 
symmetrical as possible about two mutually perpendicular planes in order 
to minimize the effect of thermal deformations resulting from changes 
in flow recovery temperature. 


Figures 4(a) and 4(b) give the calibration curves of the floating 
element under a static load. The readings are voltages induced in a 
Schaevitz coil (see fig. 5) by motion of tie core due to the loading. 


Models 

The flat plate used in the 8- by 8- inch tunnel (figs* 1 and 6), 
made of smoothly finished ground die steel, completely spanned the tun- 
nel. Rubber gaskets were glued to the edges to prevent air leakage 
around the sides and to maintain two-dimene ional flow. Forty-six 
static-pressure orifices gave chordwise anc. spanwise pressure distribu- 
tions which were recorded on a mercury manometer board. An iron- 
constantan- thermocouple junction located at the same chordwise station 
as the element center line provided wall temperature readings. 



7 


Figures 2 and 7 show the flat-plate model used in the 2- 


by 2— - inch 
J 2 


tunnel which was also made of ground die steel. The plate was firmly 
held flush against one of the tunnel walls with its sharp leading edge 
extending just beyond the boundary- layer -removal scoop. There were 
13 static-pressure taps on the surface. A row of eight holes, spaced 
l/4 inch apart, located 0.30 inch from and parallel to the leading edge 
served as a boundary- layer trip. These holes were interconnected with a 
brass tube which had an outside lead equipped with a throttle. The pres- 
sure differential between the test section and the room induced a mass 
flow sufficient to trip the boundary layer. Two-dimensionality was main- 
tained by O-rings inserted in channels machined along the sides of the 
plate. Wall temperature was measured with an iron-constanton 
thermocouple . 


Stanton Tube 


A sketch, a photograph, and a photomicrograph of the Stanton tube 
axe shown in figures 8(a), 8(b), and 8(c), respectively. It was made 
of a hypodermic tube flattened and honed so that maximum probe -tip 
height and opening were 0.004 and 0.001 inch, respectively, as deter- 
mined by measurements with a 62:1 optical comparator and a 30 -power 
microscope. Both tunnels were equipped with traversing mechanisms which 
allowed motion streamwise and normal to the plate. The difference 
between the probe pressure and the appropriate static pressure could be 
recorded on either a mercury or a di-n-butyl phthalate U-tube manometer. 
The response time of the di-n-butyl manometer was of the order of 
3 minutes. 


The distance from the probe to the flat plate was determined visu- 
ally by using a telescopic cathetometer (fig. 9). It is believed that 
the probe location could be determined well within ±0.001 inch. It was 
also found that even when the tube rested on the surface, as was required 
for the shear calibration, its lower edge did not quite touch the plate. 
The center line of the probe tip, when resting against the flat surface, 
was 0.0035 snd- 0.0028 inch away from the wall in the 8- by 8- inch and 


2- by 2— - inch tunnels, 
2 


respectively. 


EXPERIMENTAL RESULTS 


The experimental program consisted of obtaining direct shear meas- 
urements with the floating element and pressure readings with the 
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Stanton tube . The skin-friction measurements served as the reference 
values against which the probe was calibrated. 


Floating -Element Shear Measurements 

Figures 10 and 11 show the skin-friction-coefficient values measured 
in the 8- by 8- inch and the 2- by 2|- inch tunnels, respectively. The 

laminar reference curve used in figure 10 is calculated frcm Chapman and 
Rubesin's theory (ref. 14). The turbulent “cc-pressible reference plots 
are the semiempirical curves given by Von Karman (ref. 15) and Coles 
(ref 16)! It is readily seen from figure 10 that the subsonic runs 

were in transition. 

The laminar supersonic data (fig. 10) are represented rather well 
by a curve parallel to the theoretical one and 10 P e J-^nt hugher The 
standard deviation due to scatter of the c.ata about the be: st ^ 
is 11 percent as computed by the least- sqi are error method. The scatter 
is believed to be caused mainly by drift in tunnel conditions which may 
affect the Reynolds numbers. Variations in the instrument null pos 
So contributed to errors. The subject of gap effects is covered 
thoroughly by Hakkinen (ref. 2) . 

The turbulent skin-friction coefficients (fig- ll) axe plotted 
against R 0 , the Reynolds numbers based on the local momentum thickness. 
This is a more significant parameter than B* since the transition pom 
is not well defined. Previous investigat .ons have not presented turbu- 
lent 0 shear measurements in supersonic flow for values of R e below 1,500 
Ac; a matter of fact Coles (ref. l6) questions the validity of the 
turbulent-boundary-layer similarity concept for R 0 = 2,000. The presen 
data vary from Re - 6 5 0, which is very near the critical Reynolds number 
to R = 2 500 The skin-friction measurements shown in figure are 

very good Agreement with Hakkinen ’s data presented in reference 2. The 
compressibility effects do not appear to be significant below - 1.7. 

An additional source of experimental errors in turbulent shear 
measurements is the sensitivity of the mcnentum thickness 

P u L _U 

U 


0 


f (-\ _ JLVy to the velocity- profile shape near the outer 

J o P- U °° \ U °°' 


edge of the boundary layer. Typical turbulent velocity and momentum 
profiles are shown in figures 13 and 14. 
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Stanton Tube Calibration in Subsonic and 
Laminar Supersonic Flow 

5- _ ( P t “ P o) p w 


Figure 12 is a logarithmic plot of P 
m _ V w 


against 




hr 


where P T , P Q , t w , p^, and n w are the probe pressure. 


static pressure, wall shear stress, wall density, and viscosity, respec- 
tively. The data are represented rather well by the line 


5/6 

p = o.oo6t 


o 

If P and T are nondimensionalized by division through by h , where 

h = 0.01 centimeter is the probe-tip height, the dimensionless calibra- 
tion formula is found to be 


P = 2.79-r 5 / 5 (1) 

where p = P/h 2 and t = t/ii 2 . The constant of proportionality is 
of order unity which is in agreement with the theoretical prediction 
(See section entitled "Theory of the Stanton Tube.") 

As seen from figure 12 the calibration plot is an excellent repre- 
sentation of the subsonic points which exhibit a standard deviation of 
3 percent from the curve. The supersonic points are also well described 
by the drawn line, with a scatter comparable to that shown in figure 10. 
This is to be expected since t is essentially proportional to c^. 

A plot of p against t computed from the Chapman- Rube sin formula 
shows a scatter of the order of 3 percent so that the main source of scat- 
ter in the supersonic measurements is probably in the measurements of skin 
friction. 


Stanton Tube Calibration in Turbulent Supersonic Flow 

The experimental results for the case of turbulent supersonic flow 
are not so clear cut or as amenable to analysis as the laminar and sub- 
sonic data. The best fit for the data in figure 15 is given by the line 

- — 4/3 

P = 0.390T ,J 
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or, reverting to the ditnensionless form, 


P = 7.95t 4 / 3 


( 2 ) 


The scatter is considerably larger than that exhibited by the data 
in the plot of c^ against Rg (fig- ll) - This would seem to indicate 

that not all the scatter in the Stanton tube calibration plot is due to 
the random dispersion of experimental data. It is believed that because 
of the relatively narrow Mach number spectrum, 1*5 < Moo < 1-9* the devi- 
ations due to compressibility are of the same order as those induced by 
the random scatter and hence the difficulty in separating the two effects 


Velocity Profiles 

Two groups of velocity profiles were taken in the laminar super- 
sonic runs and are shown in figure l6. Those labeled group A were 
obtained at a distance of 3*0 inches from the leading edge at Reynolds 
numbers exceeding 400,000. Group B represents profiles measured 
1.1 inches from the leading edge at considerably lower Reynolds numbers. 
They are compared with the theoretical compressible boundary velocity 
profiles predicted by Chapman and Rubesin (ref. 14). The profiles of 
group A are in reasonable agreement with the predicted ones. Except 
for Tj < 1, for which the experimental point£ coincide with the theoreti- 
cal curve, the experimental data indicate an increase of about 10 to 
12 percent in rj. Bradfield (ref. 17 ), taking profiles at a Mach num- 
ber of 3 &n<i a Reynolds number of 225,000, noted the same discrepancy. 

In the present case, however, it is seen thal the measured skin-friction 
stresses, given by the slope at r\ = 0, are Ihe same as the predicted 
ones . 


The profiles of group B do not agree as well with the theoretical 
curve. These velocity profiles were taken at lower Reynolds numbers; 
uncertainties in probe location therefore gave rise to higher percentage 
errors. Another source of error was due to the fact that these veloc- 
ity profiles were measured on the wedge part of the flat plate very 
near the leading edge (1.1 inches). Korobkir (ref. IB) has shown that 
recovery factors on the wedge are higher thar on an insulated plate. 

This effect is due to heat transfer from the lower wedge surface to the 
upper one, caused by the higher recovery temperature of the compressive 
flow at the lower surface. The same effect vas also observed experi- 
mentally by Shoulberg, Hill, and Rivas (ref. 19)- Had this effect been 
compensated for, the observed values of tj would have been reduced 
somewhat. It should be noted, furthermore, that even group B profiles 
give the correct wall shear within 10 percent . 
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Velocity profiles were taken for about half the turbulent runs, 
typical profile is shown in figure 13- The slope at Y =0 was made 
to^ match the wall shearing stress as measured by the floating element. 
The momentum thickness of this profile is calculated from the plot of 


pU A _ JL\ against Y as shown in figure 13- Notice that 0 is about 

PooUoA Uoo / 

l/86, which is typical of turbulent boundary layers. 


Corrections Applied to Laminar Velocity Profiles 

In figure l6 it is observed that the points closest to the wall 
yield high values for U/U M when the data are reduced using the Bernoulli 
relation. It was decided to use the fact that the same total-head probe 
was used in measuring the profile and the skin friction. Near the wall 
the velocity profile is linear and 


U _ 1 /V \ _ K Y \ I 
Uc UoqIm-v / IPoo^oo/h^ 


(3) 


(u/u„) c 

t is obtained from the calibration curve 


This becomes the corrected value 
pressure, 


when, for a given measured 


This correction was applied to the velocity profiles of both 
groups A and B. The corrected values are indicated in figure 16 by flagged 
symbols . The success of the correction is most notable in the case 
group A where the corrected points are seen to coincide with the theoret 
cal profile. Because of the congestion of the many test points of group 
the shift due to the correction is not so clearly visi e &s 
case of group A, but all the corrected points of group B lie near the 

theoretical curve. 

It is therefore suggested that the same correction is applicable 
even in the case where the probe had not previously been calibrated 
against an absolute skin- friction meter. The procedure would be modified 

somewhat. In the relation p = Kt 5/3, K is the constant of propor- 
tionality characteristic of the probe and is determined from a zero- 
pressure-gradient velocity profile by letting t be the theoretical 
vSue With this information the "high" points of other profiles can be 
corrected. This special procedure is most useful in cases where the 
profiles to be corrected are measured in regions for which a theoretical 
shear value is not easily available (e.g., decelerating flows over air- 
foils, etc.). 
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THEORY OF THE STAN! ON TUBE 


The theory of the Stanton tube rests on the assumption t ha t its 
characteristic height h is small compared with the boundary -layer 
thickness (h/6 << l) . Within this assumption there are two classes of 
problems. If h/6 is small enough, the pressure and velocity perturba- 
tions induced by the probe are confined to a very narrow band in which 
the boundary- layer velocity profile is linear and the mean flow in this 
sublayer is essentially incompressible. The other class consists of the 
cases where h is still considerably smaller than 6 but the perturba- 
tion field is no longer confined to the linear portion of the profile. 
These two classes of problems will be referred to as case I and case II 


As long as the mean flow in the sublayer is incompressible the 
difference P T - P Q between the Stanton tube pressure and the wall static 
pressure depends only on the shearing stress r w , the density p v , the 
viscosity M-y., and the tube height h. As pointed out in references 2 
9*nd 10 dimensional analysis gives 

^ S P T - P 0 = TvF ^^" j = V( T ) (M 


where, by the H- theorem. 


^/ T v 


and 


T wPw h2 

p— are the only two 

Pw 


independent dimensionless parameters of the problem. In fact for the 
case of the linear profile, r is the Reynolds number based on the tube 

height h, and the velocity U T = at the edge of the tube. A dimen- 


sionless form of equation (4) is obtained 


(^ 2 Aw 2 ) 


so that 


tlirough multiplication 


by 



= P = tF( T ; = f(-r) 


(5) 


The desired calibration curve is then given by f(r). When the charac- 

t , 1C , Reyn ° ldS nUmb6r T 18 Sma11 ( T « the perturbation flow 
induced by the Stanton tube is a Stokes flow and F(t) is a constant. 
Taylor (ref. 5) estimated that constant to be 1.2. 
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When the free-stream Mach number increases sufficiently so that 
M h , the Mach number at the upper edge of the tube, is appreciable, the 

fluid compressibility becomes important. In that case equation (5) must 
be rewritten as 

P=f( T ,M h ) (6) 

Let M t be the probe-tip Mach number in case I (i.e., uniform 
shear flow). Then the parameter M T is related to the state of the 
fluid by the equation 


U 2 _ U T 2 _ (v VO: t2 

a 2 ?( P o/P w ) *Po 

where p Q is the dimensionless static pressure defined by 


( 7 ) 


P oPw h£ 


Case I - The Stanton Tube in Uniform Shear Flow 

An attempt is now made to derive an analytical expression for the 
calibration function f(-r) for the case in which the local Reynolds 
number t is large. The theory assumes that the Stanton tube creates 

a boundary sublayer the thickness of which is of order h^ « lj. In 

this inner sublayer the velocity profile is linear* and the density and 
temperature are presumed to be constant. The present a nal ysis for case I 
has been reported previously by Trilling and Hakkinen (ref. 20). 

The dimensionless perturbation equations for the sub lay er are 
obtained by perturbing the equation of a simple incompressible Couette 
flow in the manner of Orr and Sommerfeld. They are 

ys x + u x + v y = 0 


Px 1 rp 

+ v + = ~i v u + 

T*“ T L 



yv x + 
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In these equations all lengths are made dimensionless by division through 
by h, the probe-tip height; velocities, by d. vision through by T w h/p w , 

the mean velocity at the upper edge of the tune; density s, by division 
through by p . The pressure is as indicated in equation (5). 

The perturbation flow is essentially isentropic. Using the isentropic 
relation p/s? = Constant and the definition of M T , the equations become 



+ Ux + Vy 


= 0 


( 8 ) 


y% + 


V 




(9) 


y^x + 




(10) 


The boundary conditions on the lower edge of the sublayer cannot 
be specified clearly since the details of the flow structure near the 
wall are very complicated. One condition to be selected is v(x,0) = 0. 
The component u(x,0) must vanish on the wal!. up to the point x = -L 
where the streamline i|r = 0 leaves the wall; it vanishes on the surface 
of the tube; in the region -L < x < 0 it is a function of x of order 
unity the precise form of which changes with ~,ube geometry and is dif- 
ficult to specify. All perturbations are bounded on the outer edge of 
the sublayer. Thus, 


v(x,0) = 0 

(11a) 

u(x,0) s Uq(x) = 0(l ) 

(Hb) 

lim u,v < ® 

(He) 

y->°° 



Next Prandtl's boundary-layer argument is applied to equations (8) 
to (10). It is desirable to stretch the normal coordinate and the normal 
velocity component so that in the resulting x momentum equation all 
the variables and their derivatives are of order unity. The preceding 



15 


argument is valid as long as u Q (x) = 0(l) and r » l, so that equa- 
tions (8) to (10) really represent a boundary- layer -type flow. Let the 
stretching be as follows: 


I = 
n = 



(t » 1) (l2a) 


= u(x,y) 

= T k v (x,y) k 


P = rHp 


(12b) 


where k and n are the "stretching" constants to be determined so that 
the transformed equations are compatible with the order-of-magnitude 
requirements . 


The transformed equations are 


x - (2+k+n) r ^2^ + = o 


(15) 


T . T l“5k— l-3k— -(l+2k+n)_ _ -2k A- — \ / 

^ ^ + T V+T P I = U ^ + T (j U H + %) (14) 


qT- 2k 7, + r k - n - 


n " 2 P1 = T" (k+1 )v + T 1-k /^- v + - n 

55 pm 3 1 




(15) 


If the inertia terms on the left-hand side of equation (14) are to balance 
the viscous terms, the exponents of t must vanish in all the inertia 
terms in equation (lL); that is, 


1 - 3k = 0 


1 + 2k + n = 0 



l6 


from which 


k = 


1 

5 


n = 


5 

3 


Substituting these values of k and n ini o the equations of motion, 
equations ( 13 ) to ( 15 ), gives 


t + u_ + = 


0 


( 16 ) 


i u s + v + h - "tih + t ' 2/5 | u ii + v j n ) 


(17) 


T _ 2 / 5 qV| + p 


T ' 2/5 (l ^ + 1 %) + T ' V ^i (18) 


Since t » 1, all terms in t" 2 /^, and so forth are negligible 

in comparison with terms independent of t. It is seen from equation (l8) 


that p^ = o(t' 2 /3) « which is the usual boundary-layer result of a 

negligible normal pressure gradient across tie layer. The approximate 
sublayer equations are then 


u^ + v q = 0 


(19) 


qu t + V + p (6) = u^ (20) 

Referring again to the transformation equation (l2b) it is seen that 

- -n - 5/3 
P = pr = pT 


and since p = 0 (l), the required form of f("r) in equation ( 5 ) is 
found to be 
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p = At5/3 

where A is a constant of proportionality of the order of 1 . 


( 21 ) 


Case II - The Nonlinear Profile 

As t increases to very high values, either because of the large 
size of the Stanton tube or because of the high shear in a turbulent 
flow, it is possible that the probe, and hence the perturbation field, 
will extend beyond the linear portion of the velocity profile. In that 
case, the previous analysis which led to the 5 / 3 -power law is invalid. 

The physical situation is shown in the following sketch: 



The calibration depends strongly on the shape of the profile. To 
simplify the analysis the actual physical situation is replaced by the 
following model: 
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The constant velocity is taken to be proportional to the turbulent 

friction velocity U* = /’VAV According to this model the form of the 
calibration function depends on the ratio of sublayer thickness to tube 
height, a = e/h. 

The problem can now be broken into its component parts, linear pro- 
file region and constant profile region. The effect of each can be 
found separately and then the two are superimposed with weighting func- 
tions cl and 1 — cc being multiplied intc the linear and constant 
profile contributions, respectively, as 

p = aAf ( t ) + (l - a)Bg(-r) (22) 

where A is the constant appearing in equation (2l), g( T ) is the 
characteristic function of the calibration curve for a constant veloc- 
ity profile, and B is the constant of proportionality associated 
with g(-r). The justification of this procedure rests on the fact that 
the two terms in equation (22) add up to give the perturbation pressure 
at Y = h. 


Next g( t ) is evaluated, f(-r) being found in the preceding sec- 
tion. The perturbation equations correspording to equations (8) to (10) 
sire: 



+u^+ v- 


= 0 


(23) 


r c u x 


Px 

T 


1 _ 

fiF 





( 21 +) 


T 


e v x 



\F 



+ 



+ 



(25) 


In these equations the nondimensionalization is performed as before 
except that the velocity at the upper edge of the probe is taken to be 

t £ U* instead of U* = ^r w /p w is the friction velocity connected 

with the sublayer of a turbulent boundary layer; and t £ is a dimen- 
sionless universal constant for all turbulent boundary layers. It is the 
velocity (made dimensionless through division by U*) of the flow at 
the outer edge of the laminar sublayer or the Reynolds number of the 
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sublayer. The numerical value of t £ is of the order of 20 . Also 
M* 2 = U* 2 /a 2 = t/tp 0 . 

If the same stretching procedure as before is followed, the trans- 
formed equations are 


t £ M* 


2 T -(n+l)„ 


p 6 + 


u e + v n 


= o 


(26) 


1 

t2 


-2k 


T e u ^ 



T -2k/iu + ly ] 
\3 ^ 3 in] 


(27) 


-k — 
T T e v | 


. _k-l— 
+ T P„ 




+ 


4 

3 


T k y 


nn 



( 28 ) 


from which k = l /4 and tj = -1, so that g("r) in equation (22) is 
given by 


g(*0 = t 


and equation (22) becomes 


= oAt ^ 5 + (l - oc)Bt 


p = OAT"' + \1 - ajilT ( 29 ) 

Next it is observed that a can be expressed in the following manner: 


a = e/h = 


( u VAw) 




( 30 ) 


Substituting equation (30) into equation (29) gives 


P = At 


t 5/3 


+ Bt 


^1 - = At £ t 7 / 6 + b(Vt - T e )^ (31) 
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For large values of t (say t > t £ 2 ) the first term on the right- 
hand side of equation (31) predominates so that 


p At 


, t 7/6 „ a . t 7/6 


( 32 ) 


where A' is not of order unity but rather of order t £ . 

This relation can be obtained in a different manner which is more 
useful for an extension of the theory to the compressible case. For 
a- 4 0, it can be deduced directly from equation (29) that for sufficiently 
high values of t 


P ~ 


OAT?/? = 


apj. 


( 33 ) 


where p^. is the calibration function for ;ase I. An evaluation of a 

may be made using an argument that reconciles the actual velocity profile 
and the assumed model. Consider the geometry of the case shown in the 
following sketch: 



In this sketch is the actual velocity at the upper edge of the probe 

and U T is the velocity that would have existed at the upper edge of 
the probe if the profile had continued to b€ linear beyond e. Then, 


a = 
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so that equation ( 33 ) becomes 



where is the desired calibration formula in case II. 


( 34 ) 


It will now be demonstrated that Uj 1 / AJ t is indeed proportional 
to i/|/t\ The actual velocity at h is U^. For large values of t, 
Ujj can be evaluated from the logarithmic law for turbulent velocity 

profiles. Van Driest (ref. 21), for instance, gives for a smooth wall 
(using the present notation) 


a 

. = _ = 5-24 + 5.75 log \fF 


The velocity U T is defined as 

U T * 

so that 



^Jp v (3.24 + 5-75 log vF) 


3.24 + 3.73 log \fr_ 

\F 


( 35 ) 


2 6 

As long as r is neither too small nor too large, say 10 < t < 10 , 

the behavior of Uh/U T is dominated by the denominator of equation ( 35 ) » 
that is, 

(»h/«r) « ~ (36) 

Combining equations (34) and (35) gives 



oc 




A"t 


7/6 


( 37 ) 
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1 1 ^ 

It should be noted that A is not of order unity since the numerator 
of equation (35) is about 20 ^i.e., of order T e) # Thus, equations (37) 

and (32) agree both in form and magnitude. 


COMPARISON OF PRESENT AND PREVIOUS EXPERIMENTS WITH THEORY 


Trilling and Hakkinen (ref. 20 ) gathered and presented systematically 
a wide variety of Stanton tube calibrations in subsonic flow. It is 
found that for t < 600 the calibration functions follow the 5/3-power 
law. The present subsonic measurements fall within this range and agree 
with the previous results. The data of Preston (ref. 10) for t > 600 
and those of Fage and Sargent (ref. 6) for r > 200 deviate appreciably 
from the 5/3-power law. In Preston's investigation p is found to be 

proportional to t F age and Sargent obtained a somewhat higher loga- 
rithmic slope for their measurements. It is suggested that in both 
instances the calibration falls within the case II classification and 
that the reported results are, within experimental errors, in agree- 
ment with the predicted 7/6-power law. In the case of Preston’s data 
the difference in the logarithmic slopes is less than 3 percent - 
certainly within experimental accuracy. The constants of proportional- 
ity in the data of references 10 and 6 are 33 and- 18, respectively. 
Allowing for differences in tube geometries, these values agree well 
with the predicted constant A’ of order r e , which is about 20. 

The present results in the laminar supersonic case seem to indicate 
that for moderate Mach number and laminar flow, compressibility effects 
are unimportant. 

The turbulent supersonic results, yielding a 4/3-power law, can be 
explained by extending the results of the section entitled "Case II - 
The Nonlinear Profile." For the incompressible flow case the following 
shear correction was obtained: 


p “ p u = (^) p i <5 

and. since U^/Ur was found analytically to be proportional to l/\/T 
the 7/6-power law resulted. If the argument :.s extended by analogy to 
the compressible case, the relation is expected to be of the form 


P = 



P 


I 


(39) 
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where is the actual Mach number at the upper edge of the probe, 

2 / * 

M x is t j 7p^, the Mach number which would have been obtained if the 

profile were linear up to h, and pj is the calibration formula for 
a linear profile. 

Unfortunately, the Mach number profile for a supersonic turbulent 
boundary layer is not known analytically so that the dependence of 
Mh/^T on t is not known. The alternative is to determine ^ 1 /M t 

experimentally. Figure 17 presents a plot of M T /M h computed from 

velocity profiles and static-pressure measurements. The data are well 
represented by the line 


M T /M h = 0 . 18 t 1/5 ( 40 ) 

Substituting equation ( 40 ) into equation (39) gives 



as compared with the experimental calibration value of 

P = 7 . 95 t 4 / 5 ( 42 ) 

Apparently two compressibility effects axe encountered. Not only is 
the calibration power law changed through the Mach number dependence 
of the shear-correction factor but the proportionality constant 

is also affected directly. It is lower than expected from equation (32). 
Hakkinen (ref. 13) found at M^ = 1.5 that p/p t = p/p = 2.3. 

7 ^subsonic */ v I J 

In the present investigation for the range of Mach numbers covered 

(1.3 < m*, < 1.9) P /p T = L 22 . = 2.85. 

/ 1 2.79 

SUMMARY OF RESULTS 


Available information on the calibration of the Stanton tube as a 
skin-friction meter has been reviewed and analyzed. Additional experi- 
ments were performed in flow regions not covered in previous 
investigations . 
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A theoretical model is developed to explain the subsonic results in 

the range 10 < t < 10^ of the dimensionless shear parameter. This model 
is shown to be adequate also for laminar supersonic flows. A shear cor- 
rection applied to the basic model is vised to predict the deviations 
occurring at higher values of t when the perturbation field extends 
outside the linear portion of the boundary -layer velocity profile. 

The main results are shown to be: 

(1) In subsonic flow, for 10 < t < 10^, the dimensionless pressure 
is related to the shear parameter t by an equation of the form 

p = Kt5/5 

where K is a constant of the order of unity which depends on the tube 
geometry only 

(2) These results hold also for the case of laminar supersonic 
flow for Moo = 2 and. t = 100. The upper limits of M^, and t in 
the 1 nfivl nflr supersonic case for which this result is still valid have 
not been established. 

(3) In subsonic flow, for t > 600 and for values of h/6 such 
that the probe sticks out of the linear portion of the profile, the 
calibration formula is 


P = 



where K' is of order t £ , the universal turbulent Reynolds number 

based on the friction velocity ]j T v / Pw an ^- sublayer thickness. 

(4) In supersonic turbulent flow ( 1. 5 < M» < - 1*9) the calibration 
function is of the form 


where K" 
effects. 


is smaller than 


K 



, apparently because of compressibility 


All the results concerning the form of rhe calibration function 
just summarized may be explained by the theoretical model and the shear 
correction. 

No compressibility effects were detected in the laminar supersonic 
tests. The range of Mach numbers in the turbulent supersonic experiments 
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was insufficient to isolate the compressibility effects. Further inves- 
tigations over a larger range of Mach numbers and shear parameter t are 
needed before the compressibility effects can be fully understood. 

A study to investigate the implication of the rapid density varia- 
tions in the supersonic turbulent boundary layer is desired. It is 
possible that t, as defined herein using wall conditions only, is not 
the best -suited form of the shear parameter. 


Massachusetts Institute of Technology, 
Cambridge, Mass., July 1, 1957. 
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L- 58-3992 

Figure 1.- Flat plate mounted in 8- by 8- inch supersonic wind tunnel 
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L-58-3995 

Figure 2.- Flat* plate installation in 2- by 2-- inch supersonic wind tunnel. 
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(a) Sketch showing cutaway view. Scale, 1 inch = 1/2 inch 
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(a) Laminar experiments. 

Figure 4.- Calibration of skin- frit tion meter. 
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(b) Turbulent experiments. 
Figure 4 . - Concluded . 





Figure 5*- Electrical circuit. 
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Figure 8.- Stanton tube configuration. 
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Figure 8,- Continued. 




(c) 


^ L-58-399 6 

Photomicrograph of* head-on view of probe tip. 


Figure 8.- Concluded. 
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Figure 9 -- Installation of cathetometer . 
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Figure 12.- Stanton tube calibration. 
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Figure 15.- Stanton tube calibration in turbulent supersonic flow. 
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Figure 17 •- Deviation from linear profile in turbulent flow as a function 

of shear parameter t. 
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